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The many sources effect on tfie genuine multifiadron correlations * 

Gideon Alexander''^ Q^^id Edward K.G. Sarkisyan''^ 

'^School of Physics and Astronomy, The Raymond and Beverly Sackler Faculty of Exact Sciences, 
Tel- Aviv University, IL-69978 Tel- Aviv, Israel 

Here we report on a study aimed to explore the dependence of the genuine multiparticle correlations on the 
number of sources while the influence of other possible factors affecting the multihadron production is avoided. 
The analysis utilised the normalised cumulants, calculated in three-dimensional phase space, of the reaction 
e^e~— ^Z" — > hadrons using a large Monte Carlo event sample. The multi-sources reactions were simulated by 
overlaying a few independent single e'''e~ annihilation events. It was found that as the number of sources S 
increases, the cumulants do not change signiflcantly their structure, but those of an order q > 2 (i.e. more 
than 2 pions) decrease fast in their magnitude. This reduction can be understood in terms of combinatorial 
considerations of source mixing which dilutes the correlations by a factor of about l/S"*^^ which can also serve as 
a method to estimated the number of sources. This expected suppression is well reproduced by recent cumulant 
measurements in hadron and nucleus induced reactions both in one (rapidity) and two (rapidity vs. azimuthal 
angle) dimensions. The diminishing genuine correlations effect should also appear in other dynamical correlations 



like the Bose-Einstein in e"'"e — >W+W 



hadrons and in nucleus-nucleus reactions. 



1. INTRODUCTION 

During the last decade an increase interest has 
been shown for the genuine multiparticle corre- 
lations in multihadron final states of hadronic, 
e"'"e~ and other reactions §. Recently OPAL, m 
its study of e+e" annihilations on the Z*^ mass, 
has established the existence of strong genuine 
multihadron correlations up to the fifth order . 
In hadron-hadron, like proton-proton, collisions 
the correlations of more than three particles have 
also been observed 1^-^. In contrast to this sit- 
uation, in heavy ion collisions, at low energies 
and/or in reactions of light nuclei, genuine corre- 
lations are found to have non-zero values only up 
to the third order Furthermore it has been 
found out that in general these correlations be- 
come weaker as the reaction average multiplicity 
increases. In nucleus-nucleus collisions at high 
energies, of tens and hundreds GeV per nucleon, 
the two-particle correlations are the only one that 
survive §|-||. 



*Invited talk presented by G. Alexander at the 9th In- 
ternational Workshop on Multiparticle Production, Turin, 
Italy, June 12 - 17, 2000. 
tEmail address: alex@lepl.tau.ac.il 
■tEmail address: edward@lepl.tau.ac.il 



This correlation dependence on the average 
multiplicity is very similar to the one observed in 
the investigations of multiparticle dynamical fluc- 
tuations, i.e. variation of many particle bunches 
in restricted phase space regions Q]. In these 
studies, known as intermittency analyses, the 
observed average multiplicity dependence of the 
correlations has been proposed to be the conse- 
quence of a mixing of several independent emis- 
sion sources |^,10,ll[. As a result, the dynam- 
ical fluctuations in nucleus-nucleus collisions are 
already well accounted for by two-particle corre- 
lations ItIP^ , whereas in hadron-hadron interac- 
tions fl, |l3| and in e"'"e~ annihilations higher 
order genuine correlations do exist. 

An analogous situation seems also to exist in 
the Bose-Einstein correlations (BEG) where iden- 
tical bosons are correlated when they emerge from 
the interaction in nearby phase space. A gen- 
uine three-pion BEG has been detected in hadron- 
hadron reactions jlj] and found to be even more 
pronounced in e+e~ annihilations |l5[|. On the 
other hand the NA44 collaboration [|16[ has not 
found, in their study of sulphur-lead collisions at 
200 GeV/A, any genuine three-pion BEG as the 
three-body correlations were well reproduced in 
terms of two-particle BEG. In a recent reprint, the 
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WA98 collaboration [0 reported on their study 
of central lead-lead collisions at 158 GeV/c per 
nucleon where an evidence for a three-pion BEG 
seems to exist. However it was pointed out that, 
at least some of the BEG analyses carried out 
with the super proton synchrotron at GERN, may 
have to be reevaluated [Q. Since the intermit- 
tency phenomenon and BEG seem to be closely 
related Q , the dependence of many sources on the 
strength of the BEG cannot be excluded. The su- 
perposition of emitters may also be a reason for 
the suppression of BEG of hadrons produced from 
W-boson pairs in e+e~ annihilations at LEP2 en- 
ergies where the overlapping of hadrons affect the 



accuracy of the W mass measurements |19 

All this, as well as the obvious intrinsic inter- 
est in the genuine correlations which carry most 
of the dynamics of the hadron production pro- 
cess, points to the need of dedicated studies aimed 
to investigate the correlation dependence on the 
number of emission sources. Here we study this 
dependence by grouping several e"'"e~ — > hadrons 
events to represent a multi-emission sources of 
particles. To obtain significant results, even when 
only few sources are considered, one needs a very 
high statistics, like that which can be supplied 
by a Monte Garlo (MG) generated sample. This 
allows to minimise the calculation error and thus 
be sensitive to the correlations dependence on the 
number of sources. Another advantage of using 
MG generated events is in its possibility to gen- 
erate a multihadron sample free from contami- 
nation of other processes like, for example, e+e^ 
T^T" hadrons. 

Here we present the results of a MG study on 
the efii'ect of several emission sources on the gen- 
uine higher order multiparticle correlations |20| . 
The study was based on a generated sample of 
about 5 X 10^ events of the reaction e+e~ Z" 

hadrons which passed a full simulation of the 
OPAL detector at LEP and did reproduce rather 
well the measured genuine high order correlations 
present in the OPAL hadronic TP decay data ||]. 
Moreover, the e+c^ TP ^ hadrons annihi- 
lations should represent well the one emission 
source situation in contrast to events produced 
in hadron-hadron interactions. Our analysis on 
the dependence of many sources on the genuine 



correlations was thus carried out in a way that 
avoided effects coming from other features, like 
the multiplicity which was discussed recently in 
connection with two-particle BEG analyses in ||2^ 
and 11^. Inasmuch that final state interactions 
between hadrons coming from different sources 
can be neglected, our correlation study based on 
e+e" ^ — » hadrons annihilations may be ex- 
tended to other types of reactions since the hadro- 
nisation process is believed to rest on a common 
basis |2|. 

2. THE ANALYSIS METHOD 

The analysis is based on a generated sample 
of hadronic Z° decays using the Jetset 7.4 MG 
program [p4[ including a full simulation of the 
OPAL detector at LEP H. The MG sample also 
included initial-state radiation and effects of finite 
lifetimes. The parameters of the program were 
tuned to yield a good description of the measured 
event shape and single particle distributions . 

The selection criteria for multihadron events 
used here are identical to the ones previously 
utilised by OPAL in their recent data analysis 
of multiparticle correlations Q. In particular, se- 
lected events were required to have at least five 
charged tracks each having at least 20 measured 
points in the jet chamber where the first point 
had to be closer than 40 cm from the beam axis. 
The cosine of the polar angle of the event spheric- 
ity axis with respect to the beam direction was 
required to be less than 0.7 to ensure that the 
event lies within the volume of the detector. The 
sphericity axis was calculated by using all ac- 
cepted tracks and electromagnetic and hadronic 
calorimeter clusters. 

To simulate several emission sources we did 
overlay several e^'^e^ T^ ^ hadrons generated 
events and analysed the correlations between pi- 
ous as if they were created in a single event. The 
kinematic variables are defined within each gener- 
ated e"'"e~ event with respect to its own spheric- 
ity axis. For correlation analyses of variables like 
rapidity this procedure is equivalent to the one 
where the events are rotate to a common spheric- 
ity axis. This simulates multi-sources' events pro- 
ducing hadrons, here taken to be pious, emerg- 
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ing from a common emitter. To note is that in 
this procedure the average event multipUcity is 
directly proportional to the number of sources. 
In our analysis each generated event was used 
only once to avoid correlations between differ- 
ent sources, which for this reason required a very 
large MC event sample. 

To extract the genuine dynamical g-particle 
correlations, we used bin-averaged normalised 
factorial cumulant moments, or cumulants, first 
proposed in Ref. as a tool for the search of 
genuine multiparticle correlation. 

The Cq{yi, . . . ,yq) are the g-particle correlation 
functions given by the inclusive g-particle density 
distributions Pq{yi, . . . , j/g) in terms of cluster ex- 
pansion, e.g., 

C3{yi,y2,y3) 

= P3{yi,y2,y3) ~^pi{yi)p2{y2,y3) 

(3) 

+ '2pi{yi)pi{y2)pi{y3) ■ (2) 

Here M is the number of equal bins, having a 
width Sy, into which the event phase-space is di- 
vided and the subscript (3) denotes the number of 
permutations. For simplicity we show all formu- 
lae in 1-dimensional (e.g., rapidity) phase space. 

The feature of the Cg-functions is that they 
vanish whenever there are no genuine correla- 
tions, i.e., the correlations are due to those 
present in lower orders. The correlations ex- 
tracted are of a dynamical nature since the cu- 
mulants share with normalised factorial moments 
(the intermittency analysis tool) the property of 
statistical noise suppression. 

Here we computed the cumulants as they are 
used in experimental studies, in particular we 
used the form applied in Ref. Q| namely, 

Kq = — ^"-^^^ . (3) 

The fcg'"'' factors are the unnormalised factorial 
cumulant moments, or the Mueller moments pSf, 



calculated for the mth bin. These factors repre- 
sent the correlation functions Cq integrated over 
the bin and N,n is the number of particles in the 
mth bin summed over all the J\f events. The 
definition takes into account the non-uniform 
shape of the single-particle distribution and the 
bias when the cumulants are computed at small 
bins. 

The cumulant calculations were performed in 
the three-dimensional phase space of the kine- 
matic variables commonly utilised in this kind of 
studies namely: 

• The rapidity, y ^ In ^ {E + p\\) / (E - p\\) , 
with E and p\\ being the energy and lon- 
gitudinal momentum of the hadron within 
the interval -2.0 <y < 2.0; 

• The transverse momentum in the interval 
0.09 <PT< 2.0 GeV/c; 

• The azimuthal angle, < < 27r, cal- 
culated with respect to the eigenvector of 
the momentum tensor having the smallest 
eigenvalue, in the plane perpendicular to 
the sphericity axis. 

These variables are defined with respect to the 
sphericity axis, in a way and within the intervals 
similar to those used in a recent OPAL analysis |^] 
and in other cumulant studies jl]. 

3. GENUINE CORRELATIONS AND 
THE NUMBER OF SOURCES 

3.1. Monte Carlo studies 

In Fig. |l| we reproduce the cumulants of or- 
ders q — 2, 2> and 4 measured by OPAL and 
compared with those based on two MC models, 
Jetset 7.4 and Herwig 5.9 (|]. The cumulants 
were calculated in 1-dimensional sub-space of ra- 
pidity, in 2-dimensional rapidity vs. azimuthal 
angle sub-space and in 3-dimensional phase space 
of rapidity, azimuthal angle and transverse mo- 
mentum. The measured cumulants are seen to 
have large non-zero values thus inferring strong 
genuine multiparticle correlations down to very 
small bin sizes. As was also concluded by OPAL 
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Figure 1. Cumulants of order q — 2, 3 and 4 as a function of M^^^, where M is the number of bins 
of the Z?-dimensional sub-spaces of the phase space of rapidity (y), aziniuthal angle ($) and transverse 
momentum (pr), compared to two MC models. The data and the MC predictions are taken from Ref. 1^. 

Q, the high order fluctuations cannot be repro- quire also high-order correlations to be present, 

duced by lower-order correlations and indeed re- Thus higher order correlations do play an impor- 
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Figure 2. The MC predicted cumulants of order q = 2, 3 and 4 as a function of A/^/^, where M is the 
number of bins of the D-dimensional sub-spaces of the phase space of rapidity (y), azimuthal angle ($), 
and transverse momentum (pt)- The sohd symbols represent the cumulants for a single source, while the 
open symbols are the cumulants values of seven sources. 



tant role in the hadronic Z*^ decays processes. 

As can be seen from Fig. |l|, although the MC 
cumulants slightly underestimate the data start- 
ing at intermediate bin size, they reproduce well 
the over all behaviour of the cumulants as a func- 



tion of the bin-size. This fact is utilised here for 
the study of the genuine multiparticle correlations 
dependence on the number of sources. 

In Fig. 1^ wc compare the MC based cumu- 
lants of orders q = 2, 3 and 4 calculated from 
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Figure 3. The dependence of the averaged 1- 
dimensional rapidity cumulants K^^ of order q = 
2, 3 and 4 on the number of e+e" sources. The 
cumulants were averaged over the A/-range where 
they maintain an almost constant value (see Fig. 
1). The lines represent the expected dilution ac- 
cording to Eq. (^) where q is neglected in com- 
parison to the multiplicity n. The striped areas 
are the allowed regions when q is not neglected 
with respect to n (see text). 



a single e+e^ Z'^ — > hadrons events (solid 
symbols) with those obtained by overlaying seven 
such events to represent seven hadronic sources 
(open symbols). The calculations were performed 
in the 3-dimensional phase space of rapidity, az- 
imuthal angle and transverse momentum as well 
as in its projections onto 1-dimensional rapidity 
and onto 2-dimensional rapidity vs. azimuthal 
angle sub-spaces. 

The following observations can clearly be made 
from Fig. ||. 

• The existence of a dynamical component, 
i.e. rise of the cumulants with increasing 
number of bins Af , is seen to be present 
both in the single source as well as in the 
case of many sources. Although the slopes 
of this scaling behaviour are smaller for sev- 
eral sources than for a single source, they 
are still strongly present. It is also evi- 
dent that the scaling character is kept as 
the number of sources increases. For exam- 
ple, the single-source 1- and 2-dimensional 
cumulants level off at the same M values 
as those for seven sources. No such satu- 
ration exists for the one and seven sources 
cumulants in three dimensions. 

• The genuine dynamical correlations, mea- 
sured by the cumulants, significantly de- 
crease as the number of sources increases. 
This decrease is stronger for higher order 
correlations. Whereas the two-particle cu- 
mulants suffer a reduction of an order of 
magnitude as the number of sources in- 
creases from one to seven, the four-particle 
cumulants diminish by three or four orders 
of magnitude. 

• The hierarchy of the Kq cumulants is re- 
versed as the number of sources increases. 
The cumulants derived from the single- 
source events increase with increasing q- 
order so that k'^^ < K^^^ < if^^-* , whereas 
the hierarchy in the cumulants calculated 

(7) 

for seven sources is reversed namely, K2 > 

(7) (7) 

Kj^ ' > Kl ' . In addition, the multi-sources 
cumulants of order q > 2 have almost 
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(7) 

the same reduced value namely, Kj^ ' « 
kP < O(O.l). This last feature does not 
change as the dimension increases. 

• The overall dominant feature of the analysis 
results is the diminishing value of the higher 
order cumulants as the sources number in- 
creases leaving the K2 to be the dominant 
genuine multiparticle correlation. 

To analyse further the observed diminishing ef- 
fect, the cumulants were averaged over the M 
region where they are seen in Fig. ^ to reach 
an almost constant value, i.e. over M > 10 in 
one dimension (rapidity) and M > 100 in two 
dimensions (rapidity vs. azimuthal angle). The 
resulted M-averaged cumulants K^^ are shown in 
Figs. ^ and^ respectively for one and two dimen- 
sions. The corresponding values are also listed in 
Table ll^ for the 1-dimensional cumulants and in 
Table H for the 2-dimensional ones together with 
the OPAL measured data cumulants |^ of single 
e"'"e~ TP ^ hadrons events. 

As Table |l| shows, the 1-dimensional data cu- 
mulants do agree with those derived from the 
MC sample for g < 3. Even though the single- 
source MC based q — 4 M-averaged cumulant 
lies lower than the measured points it is still con- 
sistent within errors with the data. This however 
is not the case for the two dimensions as can be 
deduced from the values given in Table |[ Al- 
ready for g > 3 the measured and the MC based 
2-dimensional averaged cumulants disagree. This 
is due the faster increase of the data cumulants 
compared to those obtained from MC (see Fig. 

§■ 

From the features of the Af-averaged cumulants 
shown in Figs. || and ^ one can conclude with the 
following obvious points. 
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Figure 4. The same as Fig. || but for the 2- 
dimensional sub-space of rapidity and azimuthal 
angle. 



• In one, as well as in two dimensions, the 
cumulants of order q > 2 decrease fast with 
the increased number of sources. 

• In one dimension the hierarchy is re- 
versed already for two sources and the two- 
particle correlations visibly dominate over 
the higher order ones. 



• In two dimensions the hierarchy changes 
too, but at a higher number (> 2) of sources 
since additional correlations in the second 
variable (in our case the azimuthal angle) 
play also a role. The latter correlations are 
within a jet which are suppressed when all 
the phase space is projected into one direc- 
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Table 1 



The Monte Carlo M averaged 1-dimensional rapidity K^^ cumulants obtained in the e+e" Z" — » 
hadrons reaction compared to those measured by OPAL Q in single data events. The averages were 
taken over the M region where the cumulants reached a constant value (see Fig. |2|). 



No. of 




K^^ {M > 10) 




Sample 


sources 


q = 2 








0.45 ±0.01 


0.67 ±0.04 


1.36 ±0.21 


OPAL data 


1 


0.486 ± 0.002 


0.632 ±0.017 


0.950 ±0.225 


MC 


2 


0.260 ±0.001 


0.171 ±0.006 


0.147 ±0.034 


!! 


4 


0.124 ±0.001 


0.041 ±0.003 


0.016 ±0.012 


55 


7 


0.072 ±0.001 


0.013 ±0.003 


0.003 ±0.008 


)! 


13 


0.037 ±0.001 


0.004 ±0.003 


0.001 ±0.009 





Table 2 

The Monte Carlo M averaged 2-dimensional (rapidity vs. azimuthal angle) K^^ cumulants obtained in 
the e"'"e~ — > Z'^ — > hadrons reaction compared with those obtained in a recent OPAL measurements Q 
of single data events. The averages were taken over the AI region where the cumulants reached an almost 
constant value (see Fig. ||). 



No. of 




K^" (M > 100) 




Sample 


sources 


q = 2 


q = 3 


9 = 4 




1.39 ±0.04 


7.91 ±0.36 


100 ± 52 


OPAL data 


1 


1.353 ±0.005 


6.0 ±0.2 


35.0 ±6.5 


MC 


2 


0.763 ±0.002 


1.85 ±0.03 


6.5 ±0.6 




4 


0.347 ±0.001 


0.398 ±0.012 


0.6 ±0.2 


» 


7 


0.202 ±0.001 


0.129 ±0.006 


0.11 ±0.05 


51 



tion, rapidity or azimuthal angle, as is dis- 
cussed in ||]. Thus it seems that there are 
features in the second dimension, here the 
azimuthal angle, that lead to higher order , 
g > 2, correlations which do not contribute 
to the 1-dimensional rapidity correlations. 

• At higher number of sources the dominant 
role of the two-particle correlations becomes 
more pronounced. 

3.2. Correlation dilution due to source 
mixing 

Within the procedure adopted here for the sim- 
ulation of multi-source events, it is clear that if 
a genuine correlation exists it can only be de- 
tected in groups of q pions emerging from the 



very same source. In those g-group combinations 
which emerge from at least two sources, genuine 
correlations should not be present. This means 
that for Kg cumulants that are calculated over 
all possible g-pion groups, the higher the number 
of sources the more diluted will be the signal for 
genuine correlations. 

For a given g-order the genuine correlation di- 
lution factor is thus: 

where denotes the number of g-particle 

groups, e.g., pairs or triplets of pions, which 
emerge from the same source. The term 
stands for the number of all possible combina- 
tions of q-particle groups which emerge from at 
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least two sources. 

Since all sources are produced in the same re- 
action and at the same energy, they do have an 
identical average charged multiplicity. For the es- 
timation of Rq we assume that all the S sources 
have the same fixed charged multiplicity n. In 
this case one has = S (p , and the dilution 
factors at g = 2, 3 and 4 are given by 



R2 



i?3 



ri>l 



(2) S + n^ (f ) 



1 

S 



2-(2)( 



l)s 



n»2 J_ 



n»3 J_ 



(5) 



where the denominators include the number of 
all possible g-particle combinations in S sources 
of charged multiplicity n. These dilution factors 
dependence on the number of sources can also be 
derived in terms of cumulants [ p^ . 

From these Rq relationes follows that as long 
as n 3> 9 one has a general expression for the 
dilution factor namely, 



Rn 



1 



(6) 



To compare the dilution factors Rq with our 
correlation results shown in Figs. |3| and U, they 
do have to be multiplied by Kq^ which is a 
measure of the genuine g-order correlation present 
in a single sources. The solid lines shown in 
Figs. H and ^ thus represent the diluted cumu- 



lants 



K] 



av(l) 



X Ra 



.q — -i-g /s-itg. The striped areas in 
which the lines are embedded are the allowed re- 
gions when q is not neglected with respect to the 
multiplicity n. 

The agreement between the cumulant calcula- 
tions and the dilution factors predictions is re- 
ally remarkable. In the one-dimensional case the 
predictions follow the MC based cumulants for 



q = 2 and 3 and are certainly well within the 
rather large errors of the g = 4 cumulants. In two 
dimensions, in addition to the remarkably good 
agreement with the MC cumulants, the dilution 
factor predictions follow closely also the hierarchy 
change of the Kq moments. 

For the order q = 2 one can relax the fixed 
charged multiplicity assumption and allow them 
to be different and still retain the Rq ~ l/S"^^^ 
relation as long as the multiplicity distribution is 
of a Poisson nature. This however is not the case 
for orders higher than 2. Nevertheless for order 
q = 3 the relation i?3 ~ derived from the 

fixed multiplicity assumption, is still valid as it 
describes well the K^^ values up to at least thir- 
teen sources (see Fig. The large cumulants' 
errors associated with the g = 4 order prohibits 
to judge how accurate is the R4 — relation. 

An additional interesting and useful applica- 
tion of the relation Rq ~ 1/S''~^ is that it offer a 
method to estimate the average number of sources 
(S) via the cumulant averaged values over the 
large M region of two sequential g-orders through 
the ratio, 



3.3. 



.(1) 
9+1 



v(l) 



(7) 



Comparison with hadron and nucleus 
induced reactions 



As is already mentioned in the introduction, 
the genuine correlations measured in e+e~ anni- 
hilations [H are found to be weaker in hadronic 
interactions 10-0, |l^ and even more so in nuclear 
collisions [|[||-||. In nucleus-nucleus collisions at 
ultra-relativistic energies only the second-order 
correlations are found to have non-zero values in 
rapidity , while in two dimensions (rapidity 

and azimuthal angle) the third-order cumulants 
arc also detected |^,^. 

In Table ^ we list the results obtained by sev- 
eral experiments on the M- averaged ra- 
pidity cumulant values for q — 2 and 3. In Ta- 
ble |4| the analogous cumulants for g = 2, 3 and 
4 are given as measured in two dimensions, ra- 
pidity vs. azimuthal angle |^,^. These average 
values were taken over the M-region where the 
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Table 3 

The M averaged 1-dimensional rapidity cumulants K^^ of orders q — 2 and 3 measured in several 
hadronic reactions. The cumulants were averaged over the M regions where they were seen to approach 
a constant value. The quoted values were estimated from the relevant published figures in the papers 
listed in the table references. 



("-ch) 






Reaction 


q = 2 


9 = 3 


~ 8 


0.32 ±0.02 


0.26 ±0.12 


7rp 


21.1 


0.21 ±0.04 


0.14 ±0.18 


pEm 


> 50 


0.34 ±0.02 


0.12 ±0.03 


AuEm 


73.3 


0.21 ±0.03 


0.05 ±0.07 


SiEm 


81.1 


0.20 ±0.05 


0.00 ±0.12 


OEm 


154.9 


0.11 ±0.05 


0.01 ±0.18 


OEm 


216.1 


0.25 ±0.05 


0.08 ±0.10 


SEm 


272.6 


0.09 ±0.05 


0.00 ±0.18 


SEm 


289.8 


0.11 ±0.08 


0.02 ±0.10 


SAu 


355.0 


0.08 ±0.05 


0.00 ±0.08 


SAu 


383.9 


0.07 ±0.04 


0.00 ±0.05 


SAu 



Beam energy 
(GeV) 



Ref. 



i3r~ 

131* 
161 

16,|] 

191 

131* 

191 
131* 

1^1** 

131* 



250 

200 

10.6A 

U.5A 

60A 

200A 

200A 

200A 

200A 

200A 

200A 



* Calculations are based on the measured factorial moments, 
collisions. 



Semi-central collisions. 



Central 



Table 4 

The M averaged 2-dimensional (rapidity vs. azimuthal angle) cumulants K^^ of orders q = 2, 3 and 
4 measured in nucleus-nucleus collisions. The cumulants were averaged over the M regions where they 
were seen to approach a constant value. These quoted values were estimated from the relevant published 
figures given in the references listed in the table. 



{rich) 



1^. 



Reaction 



Beam energy 
(GeV) 



Ref. 



II 

1^1 



> 50 
73.3 
119.3 
216.1 



0.79 ±0.02 
0.20 ±0.05 
0.57 ±0.03 
0.25 ±0.03 



1.3 ±0.1 
0.09 ±0.03 

0.3 ±0.1 
0.07 ±0.04 



3.2 ±2.8 
0.4 ±0.3 
0.4 ±0.7 
0.1 ±0.2 



AuEm 
SiEm 
OEm 
SEm 



10.6A 
14.5A 
200A 
200A 



cumulants are seen in the published figures to 
reach a constant level. The reactions and their 
cumulants values are ordered according to their 
reported mean charged multiplicity, from the low- 
est value to the highest one. 

Tables ^ and ^ show that in hadron including 
nucleus induced reactions the two-particle corre- 
lations decrease rather fast as the mean multi- 
plicity increases both in one and two dimensions. 
However the three-particle correlations are found 
to be essentially non-existing in one dimension 
even at moderately small mean multiplicity, while 
in two dimensions these correlations are seen still 
to be non-zero. Moreover, the higher the order 



of the 2-dimensional correlation the larger the 
value is. These measurements, particularly for 
2-dimensional cumulants are in an amazing agree- 
ment with our findings, see Table | and Fig. g 

Notwithstanding the possibility that produc- 
tion of hadrons in e^e^ annihilation may well be 
simpler than in hadron induced reactions, it may 
nevertheless be instructive to relate our results to 
the measured correlation data listed in Tables ^ 
and IJ. Inasmuch that the mean multiplicity in- 
creases with the number of sources, the measure- 
ments consistent with our findings. As well as in 
nucleus induced reactions, we found the decrease 
in the two-particle rapidity correlations and the 
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absence of three-pion rapidity correlations. In 
two dimensions our results show a surviving of 
higher order correlations and a change of the hi- 
erarchy not starting from two sources, but later 
at larger number of them, the effects which one 
can see in the data in nucleus-nucleus reactions. 
All this demonstrates the effect of dilution of the 
correlations with increased number of sources. 

A quantitative comparison between our find- 
ings and the correlations in nucleus-nucleus and 
hadron induced reactions is hard to make mainly 
because of the lack of information on the values 
of iiTq^^^^ which should be derived from nucleon- 
nucleon collisions data. In addition in many of 
the existing cumulant data of the nucleus induced 
reactions have still too large errors. 

Recently the two-pion BEC have been stud- 
ied 1^ in pp reaction at centre of mass energy of 
630 GeV as a function of multiplicity by using the 
normalised cumulants method similar to the one 
used here. In that analysis it has been found that 
the correlations of the cumulants of the like-sign 
pions as well as the opposite-sign pions decrease 
with the multiplicity rich- From our analysis we 
expect the pair correlation to decrease as 1/5*, 
where S is the number of sources. This indicates 
that indeed the multiplicity is at least partially 
proportional to S. The BEC dependence has also 
been investigated in the framework of the totally 
coherent emission picture and in the quantum 
optical approach |^ where the conclusions were 
that these correlations are weaker as the multi- 
plicity increases. 

4. SUMMARY AND CONCLUSIONS 

In the present work we investigated the effect 
of many emission sources on the genuine corre- 
lations in mutihadron final state. To this end 
we adopted a procedure which should minimise 
the confusion introduced by other variables like 
charged multiplicity. For the genuine correlations 
measurement we utilised the normalised cumu- 
lant method. To simulate the situation of many 
sources event we did overlay Monte Carlo gen- 
erated hadronic events treating them as one 
event. This Jetset 7.4 MC sample of some five 
million events, tuned to the OPAL data taken at 



LEPl on the Z° mass, has previously described 
rather well the measured correlations in the e"'"e~ 
^ Z*^ — > hadrons data. We studied the cumu- 
lants in one, rapidity, dimension and in two di- 
mensions of rapidity and azimuthal angle. 

The results shown here demonstrate that the 
cumulants, obtained from a single-source events 
and from events of many sources, almost do not 
change their basic structure with the decrease of 
the width of phase-space bins. This means that 
the scaling is preserved although larger slopes are 
seen to be in the case of one source compared with 
those for several sources. 

Due to source mixing the higher order cumu- 
lants are suppressed both the 1-dimensional (ra- 
pidity) and the 2-dimensional (rapidity x az- 
imuthal angle) and diminish to zero as the num- 
ber of sources increases. In both case as S in- 
creases the hierarchy is reversed and the cumu- 
lant of the lower order q — 2 dominates. This 
happens in the one dimension case already when 
S — 2 whereas in the two dimension occurs only 
at a higher S values of around S — 5. 

The correlations observed are very well repro- 
duced by assuming that genuine correlations of 
the order q can only exist when all the q hadrons 
are emerging from the same source. Therefore 
the dilution of the genuine correlation signal is 
proportional to the ratio of the probability that 
the q hadrons will emerge from the very same 
source. From simple combinatorial considera- 
tions this probability can then be approximated 
by l/S*^^^. Thus a measurement of the correla- 
tions of two sequential orders in q can be used to 
estimate the average number of sources. 

The genuine correlations measured in hadron 
and nucleus induced reactions do follow quali- 
tatively the findings of our work. In particular 
in nucleus-nucleus reactions, where many sources 
are expected to contribute to the final hadronic 
state, the rapidity cumulants of the q > 2 or- 
ders are very small and indeed consistent with 
zero. The q = 2 order of the 1-dimensional cu- 
mulants still survive but they also are getting 
smaller as the atomic number of the nuclei in- 
creases. In two dimensions, the cumulants are 
still different from zero even at the q = A order, 
however their value decreases and a change in the 
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hierarchy takes place as the muhiphcity increases. 
The general observation that cumulants decrease 
as the multiplicity increases reassures the com- 
mon notion that the higher the multiplicity the 
larger the number of sources. 

Our results may also be useful for the under- 
standing of other types of measured correlations 
like that obtained from the Bose-Einstein inter- 
ferometry of two and more identical bosons. It 
has been previously pointed out |3^] that in the 
absence of final state interactions the BEG of the 
e+e^ W+W^ hadrons will be half of that 
of the Z° decay to hadrons. From our study 
it follows directly that the two-particle BEG, or 
any other correlations, in the two-source reac- 



tion e+e 



hadrons should be re- 



duced by a factor two as compared to that of the 
hadrons emerging from one W-boson. 
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